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We investigate non-linear, spherically symmetric solutions to the coupled system of a quintessence 
field and Einstein gravity. In the presence of a scalar potential, we find regular solutions that to 
an outside observer very closely resemble Schwarzschild black holes. However, these cosmon lumps 
have neither a horizon nor a central singularity. A stability analysis reveals that our static solutions 
are dynamically unstable. It remains an open question whether analogous stable solutions exist. 



I. INTRODUCTION 

Black holes are fascinating and intriguing objects, both 
by the presence of a horizon at the Schwarzschild radius Rs 
and by the existence of a central singularity. They arise as 
solutions of pure gravity, and one may wonder if the char- 
acteristic black-hole properties can be influenced by the 
presence of other fields. In fact, adding a massless, free 
scalar field to gravity, new spherically symmetric solutions 
to the field equations have been found, for which the hori- 
zon is smoothed out 

HHH (see also 0]). In the exterior 
region (for r > Rs) these solutions can be made arbitrar- 
ily close to the Schwarzschild solution; however, the event 
horizon is replaced by a high red-shift horizon, at which 
all metric components remain finite. It was proposed that 
these "approximate black holes" might be used in numeri- 
cal simulations of black holes as a way of regularizing the 
the horizon [5| . Unfortunately, for free scalar fields the so- 
lutions exhibit a naked central singularity and were found 
to be unstable @, 0, Hi- In this Letter, we present spheri- 
cally symmetric solutions for an interacting scalar field that 
are completely regular. For a suitable scalar potential, they 
look from the outside very similar to a Schwarzschild black 
hole, but neither an event horizon nor a central singular- 
ity exist. Observational distinction from black holes seems 
extremely difficult for many such solutions. 

In order to influence the properties of the horizon, the ad- 
ditional field must be very light, with typical mass smaller 
than the inverse Schwarzschild radius. The proposal of a 
solution to the cosmological constant problem by a time- 
varying dark energy involves a cosmological scalar field - 
the cosmon, or quintessence field [§, [l(J [HI EH (see also 
[IH). The mass of the cosmon decreases with time; for 
cosmological solutions it turns out to be of the order of the 
Hubble scale [ijj • This makes this field essentially massless 
as far as black-hole physics is concerned. The cosmon is 
thus a well-motivated candidate which could influence the 



physics of black holes. 

Cosmon lumps are bound objects, for which the scalar 
and gravitational fields combine to form non-linear solu- 
tions of their field equations. Such solutions have been 
discussed in detail in [J] for situations where the influence 
of the cosmon potential is negligible. The salient features 
are: (i) the lumps are characterized by their mass and a 
suitably defined "cosmon charge" ; (ii) the black-hole event 
horizon disappears for all solutions with non-zero cosmon 
charge; (iii) for small cosmon charge, the solutions closely 
resemble black-hole solutions when seen from the outside, 
and the horizon is replaced by a sharp transition to a region 
of high red-shift. It has been speculated that such lumps, 
perhaps also in presence of an additional matter compo- 
nent, could be regarded as a candidate for dark matter, 
if sufficiently small objects are produced sufficiently early 
in cosmology and with sufficient abundance [14j. Indeed, 
interesting solutions for a cosmon coupled to dark matter 
[h| or to neutrinos [n| have been obtained. 

In this Letter we address another issue, namely the sin- 
gularity at the center of a black hole. For this reason we 
include the effects of the cosmon potential, which adds fur- 
ther non-linearities to the system of field equations. We 
find static, spherically symmetric solutions that look from 
the outside very similar to a black hole, but are completely 
regular. 



II. FIELD EQUATIONS 

We consider the system of a scalar cosmon (or 
quintessence) field (p coupled to gravity. The action reads 



S = 



M 2 R g^ v 

— 5 df^tp d u ip - V(tp) 



(1) 
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where M = 1/v8ttG is the reduced Planck mass. In or- 
der to investigate static, spherically symmetric solutions we 
employ the metric ds 2 = —B{r) dt 2 + A(r) dr 2 + r 2 (d9 2 + 
sin 2 9 d(j) 2 ). The field equations R^ v - \g^ u R = T^/M 2 
will be solved numerically for a given scalar potential V(ip). 
The energy-momentum tensor is diagonal, with energy den- 
sity p = -T° and pressures p r = T\, Pn = T% = T^^ 
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given by 
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We shall refer to these static solutions as "cosmon lumps" . 

The ten Einstein equations reduce to a set of three or- 
dinary differential equations for the functions A(r), B(r), 
and <p(r). Those equations are supplemented by the Klein- 
Gordon equation for the scalar field, Dip = dv J- v ) , 
plicitly 



or ex- 
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These four equations are redundant, and it suffices to use 
two linear combinations of the Einstein equations in addi- 
tion to ©: 
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where mo and q$ are possible contributions localized at the 
origin. Note that for a free scalar field, i.e., in the absence 
V{ip) . (2) of a potential, the quantities m and q are independent of 
r. The corresponding solutions are singular at the origin. 
Specifically, the scalar field <p and the Ricci scalar (along 
with other curvature invariants) diverge for r — > 0. For the 
regular solutions studied in this Letter, on the other hand, 
the constants mo and qo vanish, and the mass and cosmon 
charge are generated in the region where the potential is 
non-zero. 

A minimally coupled, free scalar field is one of the sim- 
plest generalizations of a Schwarzschild black hole. The ex- 
act solutions to the Einstein-Klein-Gordon equations with 
V = were discussed a long time ago in [3, Q (similar solu- 
tions in the context of Brans-Dicke theory were studied in 
the original paper [2j ) . They exhibit a naked singularity at 
the origin. The solutions are completely specified in terms 
of m, q and the cosmological value ip^ of the scalar field 
far away from any sources. Following we introduce a 
(4) new radial coordinate x by 



dip 



(3) 



i+s 



(9) 



III. GRAY HOLES 

Time-translation symmetry gives rise to a conserved cur- 
rent J p = K\R Xfl built using the time-like Killing vector 
K x = (1,0,0,0). The corresponding conserved charge can 
be identified with the total gravitating mass of the system, 



2M 2 / cPxy/gin^ = AttM 2 lim 



2 B'(r) 



r^oo ^JAiyjBjr) 

(5) 

Here g% is the determinant of the induced three-dimensional 
metric, and = (— \[B, 0,0, 0) is a time-like normal unit 
vector. Here and below, a prime denotes a derivative with 
respect to r. Generalizing this result, the enclosed gravi- 
tating mass inside a sphere with radius r can be defined 
as 

r 2 B'(r) 



m(r) = AttM z 



VMrjBir) 



(6) 



The strength of the cosmon field can be measured by the 
enclosed "cosmon charge" 



q(r) = 4ttMV 



IB(r) 
Mr) 



(7) 



where we have introduced the dimensionless scalar field 
ip(r) — 30- ip(r). Note that this charge is not conserved in 
general, since it is not derived from a global symmetry. By 
taking appropriate combinations of the field equations, the 
mass and cosmon charge can be expressed in terms of the 
potential as 



m(r) = mo — 8n / dss 2 y/A(s)B(s)V{ip(s)), 

r dv ^ 

q{r)=q + 8n ds s 2 y/A(s)B(s) —{0{s)) , 
Jo dl P 



and write the analytical solutions in the form 

4x(l + x) 



A{X) 



B(X) 



(2 X + 1 - S) 2 ' 

X 



(10) 



1 + X / 
<£(x) = Voo +7 InS(x) • 



A useful relation is dr/dx = Rh/^AB. We define the 
"wo uld-be ho rizon" R$ — 2Gm and a related scale Rh = 
2Gy / m 2 + q 2 , both of which are determined by the mass 
and cosmon charge of the solution. The coefficients 7 and 
5 are given as 



7 = 



Rs_ 

Rh 



\Jm 2 + 



(ii) 



The special case S = 1 (or 7 = 0) leads to the Schwarzschild 
solution. All other solutions with 7 / have no event 
horizon. The black hole has become "gray" . 

At large distances (r 3> Rs), the solutions for the metric 
functions obey 



A(r) = 
B(r) = 



1 - 



Rs 



R 2 g 

R% 

V24r 3 



(12) 



Large grey holes are essentially indistinguishable from 
black holes for an outside observer. This may be demon- 
strated by assuming that the massive black hole at the 
center of our galaxy with m = (3.7 ± 0.4) • 10 6 solar masses 
is a cosmon lump with 7 = q/m of order unity. The 
star S0-16 approached the center of the black hole up to a 
distance of about 600i?s [13 • l n this case, the modifica- 
tions of the metric coefficients A and B with respect to the 
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FIG. 1: The cosmon potential (I14|l for r] = 5, K = 1, and e = 0.5 
(dark), 0.3 (medium), 0.1 (light). 



Schwarzschild metric are about 7 • 10~ 7 7 2 and 2 • 10~ 10 7 2 , 
respectively. Obviously, only very large values of 7 would 
give rise to noticeable modifications of the metric at large 
distances. 



IV. REGULAR SOLUTIONS 

The solutions for a free scalar field explored in the previ- 
ous section exhibit a naked singularity at r = 0. However, 
regular solutions exist if we allow for a non-zero potential 
for the scalar field. For regular solutions we use Taylor ex- 
pansions of the form A(r) = X^^Lo A n r n for all fields, and 
find after a straightforward analysis of the field equations 



A(r) = 1 + 



V(0o) 



3Af 2 



B(r) =B (l 



0{r) = (po + 



3M 2 
1 dV(ipo) 



(13) 
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Note that the linear coefficients Ai, B\, pi vanish, while 
Aq = 1 is fixed by the demand for the absence of a conical 
singularity. Since the field equations only depend on the 
ratio B' / B, the value of Bq remains undetermined. It must 
be set in such a way that the metric is asymptotically flat, 
B(r) — > 1 for r — * 00. Thus, for a given potential, the value 
of ipo is the only free boundary condition. 

An important observation following from the above rela- 
tions is that, as r increases away from the origin, the scalar 
field climbs up the potential, unless dV (<fo) / dip — 0, which 
would require a fine-tuning. Since we should require that 
for r — > 00 the cosmon is driven to a region with zero (or 
very small positive) potential, there are essentially only two 
types of potentials that could give rise to regular solutions: 
"slope" potentials, where the scalar potential starts out at 
some negative value at the center of the cosmon lump and 
monotonically increases toward zero (or a very small posi- 
tive value) for large r, or "bump" potentials, where between 



an initial negative value and an asymptotic cosmological 
value the cosmon field climbs over a local maximum of the 
potential. From ([5]) it follows that slope potentials always 
give rise to positive mass and cosmon charge. The same 
is true for bump potentials provided the potential is suffi- 
ciently negative inside the cosmon lump. For phenomeno- 
logical purposes bump potentials appear more attractive, 
as they provide a potential barrier against collapse of the 
scalar field into the true ground state (with negative po- 
tential). Also, in this case it is possible to achieve that the 
potential asymptotically approaches the exponential form 
V((p) oc e~ 2K ^ , which is a particularly natural potential for 
the cosmon |18| . 

A simple model for a cosmon potential with a bump is 



V(0) = -V 



. 1 — tanh nip 1 — tanh Kip 
+ S e S 



(14) 



with positive parameters Vq (typically of order M 4 ), e«l, 
and k < r\. It starts at — Vb for ip <C and approaches zero 
from above for <p 3> 0, i.e. V(ip) s» Vo e e~ 2fi *, such that this 
model may explain the dark-energy content of the universe 
on large scales. The limit e — > leads to a simple slope 
potential. Figure [1] shows the model potential for some 
parameter choices. 

Our goal is to find regular cosmon-lump solutions, for 
which at r — the scalar potential takes a negative value 
somewhere to the right of the minimum of the potential, 
while asymptotically (for r — > 00) it approaches a small 
positive value to the right of the maximum. The radial de- 
pendence of the potential within a cosmon lump introduces 
another length scale Rc besides the Schwarzschild radius. 
It is the scale beyond which the effects of the potential can 
be neglected, so that the mass and cosmon charge com- 
puted from ||5J) remain constant as r is further increased. 
For r > Rc the solutions thus converge toward the so- 
lutions (fl"2")) obtained for a free scalar field. In practice, 
we solve the field equations numerically for a chosen value 
of ipo and compute the parameter 7 at a sufficiently large 
value of the radial coordinate (r > Rc)- We then match 
onto the free solution (|12p and determine the coefficient Bo, 
which fixes the absolute values m and q of the solution. 

As a particular example we discuss the solution obtained 
with the potential parameters rj = 5, K = 3, e = 0.1 and the 
starting value ipo = —0.5 for the scalar field at the center 
of the cosmon lump. The corresponding asymptotic values 
of the mass and cosmon charge are m ~ 165 and q ~ 331 
in units of M 3 /y/Vo, which implies 7 w 1.79. Figure [2] 
shows the radial profile of the scalar field and of the poten- 
tial as functions of r in units of the Schwarzschild radius 
Rs = 2Gm. Note that the potential becomes negligible 
for r > Rc with Rc ~ 0.55i?s- The dashed line shows 
for comparison the solution obtained for a free scalar field, 
which is singular at r — 0. The presence of the potential 
removes the naked singularity and gives rise to solutions 
that are regular everywhere. The corresponding results for 
the metric functions A and B are shown in Figure [3] Once 
again the dashed lines correspond to the solutions without 
the potential. They are indistinguishable from the exact 
ones for r/Rs > 0.55. 
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FIG. 2: The scalar field ip (dark) and the potential V/V (light) 
as functions of r. The dashed lines show the scalar field obtained 
with V = 0, which diverges for r — ► 0. 




FIG. 3: Metric coefficients A (dark) and 3B (light) as func- 
tions of r. Solid lines show the regular cosmon-lump solutions 
obtained with the potential (|14[) . while dashed lines show the 
singular solutions obtained with V = 0. 



At large distances (r > few Rs) the metric closely resem- 
bles the Schwarzschild solution, but it is strongly modified 
around and inside the Schwarzschild scale Rs- We find 
no horizon, i.e. no coordinate singularity, and no physical 
singularity at the center. Note, in particular, that B(r) 
converges to a non-zero value Bq ss 0.0044 for r — * 0. In 
the Schwarzschild metric, the gravitational redshift of a 
resting source 



z(r) 



1 



B(r) 



(15) 



diverges at Rs- Our solutions feature no redshift horizon, 
but a large redshift value at the center of the cosmon lump, 
see Figure |4j This behavior might simulate an event hori- 
zon, if an observer does not measure precisely enough. By 
choice of the initial parameters the redshift at the origin 
can be made arbitrarily large, but it is always finite. 



FIG. 4: The gravitational redshift z(r) for the cases of a free 
(V — 0, dashed) and interacting (V 7^ 0, solid) scalar field. In 
the first case, the redshift diverges at the singularity. In the 
latter case, the redshift is finite everywhere but becomes large 
at the center. 



The values of the mass m and cosmon charge q of the 
solution depend strongly on the initial value ipo and on the 
overall scale Vq of the potential. In this way, a wide range 
of mass scales can be obtained in our model. Interestingly, 
the value of 7 asymptotically reaches V3 in the large mass 
limit, m > 10 3 M 3 /\/K)- 



V. STABILITY ANALYSIS 

In order to estimate the stability of the regular solutions, 
we introduce radial perturbations of the metric and the 
scalar field, which are allowed to evolve in time. The full 
time-dependent field equations are then linearized about 
the static solutions (denoted by a subscript s), using the 
ansatz 



0(t,r) = tp a (r) + 5tp(t,r) , 
A(t,r) = A s (r) +6A(t i r), 
B[t,r)=B s (r)+SB(t,r). 



(16) 



Only terms of first order in the perturbations and their 
derivatives are taken into account. In contrast to the three 
static equations, one finds four independent equations for 
the time dependence of the fields, since i?oi & n d Tq\ — tpip' 
(the radial energy-flux density) do not vanish. Following 
[6j, the four linearized field equations can be combined to 
obtain one single perturbation equation for Sip. We look 
for solutions of this equation of the form 



6<p(t, r) = e 



(17) 



If a is real, the perturbation is a small oscillation about the 
static solution. However, imaginary tr leads to a dynamical 
instability. 

Introducing a new radial coordinate g(r) obeying g(0) = 
and dg/dr = \J A s j B Sl the perturbation equation for 6<p 
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FIG. 5: The pulsation potential W(r) (dark, in units of Rg 2 ) 
features a distinctive minimum in the vicinity of Rc, which sup- 
ports a shallow bound state. The corresponding eigenfunction 
9o(r)/r (light, in arbitrary units) determines the radial profile 
of the unstable perturbation to the cosmon lump. 



implies the "pulsation equation" 
d 2 * 



do 2 



with the "pulsation potential" 
W 



(18) 
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The asymptotic behavior of this function is 
~d 2 V(<p ) V(0 O ) 



(19) 
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+ 0(r 2 ), 
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(20) 



The dark curve in Figure [5] shows the pulsation potential 
for the parameter choices considered in the previous sec- 
tion. 

Equation Q18p resembles the one-dimensional 
Schrodinger equation for a particle in the potential 
W. The eigenvalues a 2 can be calculated with the 



boundary conditions that \&(0) = 0, so that Stp remains 
regular at the origin, and ^(oo) = 0, such that the 
wave function can be normalized. The knowledge of the 
lowest eigenvalue Cq is sufficient for our purposes, since it 
determines the time evolution of the perturbations. The 
cosmon-lump solution is instable if a bound state exists, 
i.e., a solution to the pulsation equation with negative 
eigenvalue Oq < 0. For our particular model, we find 
that the negative spike of the potential, centered near 
Rc « 0.55i?s, supports a very shallow bound state, whose 
eigenvalue cr 2 m — 0.197i?^ 2 is much smaller in magnitude 
than the minimum of the potential (W m j n w — 167i?g 2 ). 
The radial profile ^&o(r)/r corresponding to this eigenvalue 
is shown by the light curve in the figure. 

We have explored a large set of parameter choices for 
the cosmon potential (|14p as well as different initial values 
po for the scalar field at the center of the cosmon lump, 
and in all cases a negative eigenvalue cr 2 . exists. We thus 
conclude that the regular, static solutions to the Einstein- 
Klein-Gordon equations obtained in this Letter using a 
non-trivial scalar potential, which a priori are not in con- 
flict with the cosmic censorship conjecture (l9| . are dynam- 
ically unstable. In this respect they share the fate of the 
known singular solutions explored in the literature @, 0, HI • 



VI. CONCLUSIONS 

We have shown that a rather simple system of a self- 
interacting scalar field minimally coupled to gravity allows 
for static, spherically symmetric solutions to the Einstein- 
Klein-Gordon equations that for an outside observer look 
very close to the black-hole solutions in pure gravity, but 
have neither a horizon nor a central (naked) singularity. 
Such regular solutions are of high conceptual interest. For 
the choice of the scalar potential considered in this Letter, 
the static solutions turn out to be dynamically unstable. 
Unfortunately, they can therefore not describe astrophysi- 
cal objects such as the black holes in the center of galaxies 
or dark-matter lumps. 

Our analysis leaves open the questions whether similar 
solutions exist, which are both regular and stable. Perhaps 
such solutions can be found by modifying the cosmon po- 
tential, or by adding a matter component coupled to the 
cosmon field. If such regular "gray holes" without event 
horizons exist, this would revolutionize the conventional 
view of black holes and related puzzles such as the infor- 
mation paradox. 
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